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Abstract. The (g, i;)-Fock space ,J?'5,t(,i^^), introduced in this paper, is a deformation of the q-Fock 
space of Bozejko and Speicher. The corresponding creation and annihilation operators now satisfy 
the commutation relation 

aq,t{f)ag,t{gy ~ qag^tigYaq^tif) = {f,g)^^ , 

with N denoting the usual number operator, and generate a Hilbert space representation of the 
Chakrabarti- Jagannathan deformed quantum oscillator algebra. The moments of the deformed field 
operator Sq^t{h) := aq^t{h) + aq^t{h)* , the present analogue of the Gaussian random variable, are 
encoded by the joint statistics of crossings and nestings in pair partitions. The restriction of the 
vacuum expectation state to the (g, f)-Gaussian algebra is not tracial for t ^ 1. 

The q = < t specialization yields a new single-parameter deformation of the full Boltzmann 
Fock space of free probability. This refinement is particularly natural as the probability measure as- 
sociated with the deformed semicircular element turns out to be encoded via the Rogers-Ramanujan 
continued fraction, the t-Airy function, the t-Catalan numbers of Carlitz-Riordan, and the first- 
order statistics of the reduced Wigner process. 
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1. Introduction 

Non-commutative probability refers, in broad terms, to a generalization of the classical proba- 
bility theory to encompass quantum observables. In particular, Kolmogorov's probability triple 
is replaced by a non-commutative probability space {£/,ip), where £/ is a unital *-algebra whose 
elements are interpreted as non-commutative random variables and the unital linear functional (/? 
plays the role of the expectation. Particularly relevant examples are to be found among algebras of 
linear operators on the Bosonic and Fermionic Fock spaces together with the corresponding vacuum 
expectation states. The Bosonic and Fermionic annihilation operators, {af} and {a~} respectively, 
together with their adjoints (the creation operators) satisfy the commutation relations 

{at){ap* - (alnat) = (CCR) 

{a;){ajr + {ajnar) = d,^j, (CAR) 

where (CCR) is referred to as the canonical commutation relation^ and (CAR) as the canonical 
anti- commutation relations. 
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measure. All other examples throughout this paper will take place in *-algebras of bounded operators. 
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In the setting of (CCR), of particular note is the fact that the Bosonic field operator si = ai + al 
can be naturally identified with the classical Gaussian random variable It is in this context that 
Frisch and Bourret [1] first considered a deformation of the canonical (anti-)commutation relations, 
as a^a* — qa*ai = jijSij, for q £ [—1, 1] and ^ij some positive-definite function, and studied the 
properties of the "parastochastic" random variable aj + a*. However, it was only two decades 
later, in Bozejko and Speicher's independent study of deformed commutation relations that 
the existence question was resolved and such processes constructed. (For existence proofs, see also 
[HllllISllllITlIHlin].) In particular, Bozejko and Speicher constructed a suitably deformed Fock space 
on which the creation and annihilation operators (analogues of those on the classical Fock spaces) 
satisfied the deformed commutation relation: 

aia* - qa*ai = 6ij. 

The setting for these operators is provided by a direct sum of single particle spaces 

^ = (ra)e0jrj?'^ 

n>l 

with Q being some distinguished unit vector and the complexification of some real separable 
Hilbert space J^. On this space, one can define the deformed sesquilinear form (^,r/)g = {£,,Pqr])o, 
where ( , )o is the usual inner product on the full Fock space (cf. Section [s]). The particularly 
insightful choice of a positive "projection" operator Pg, which would turn { , )g into a bona fide 
inner product and have the classical creation operator and its resulting adjoint satisfy (g-CR), 
was the key element in Bozejko and Speicher's construction. Consider the unitary representation 

vr I— )• uir^ of the symmetric group Sn on given by ui^^ hi . . . (>S) hn = h^^ii) (8) . . . (X" h^{n)- 

Writing Pq = 0^=0^9"^ with Pg^"^ : ^ JT®", the key idea of Bozejko and Speicher was to 

replace the (anti-)symmetrization step that would yield the two classical inner products by instead 
"g-symmetrizing" as 

:= ^ (2) 

where inv(7r) gives the number of inversions in a permutation vr (cf. Section [1]). By showing that 
Pq was indeed positive definite for |g| < 1 (or positive semi-definite for \q\ < 1) and completing 
(and separating) with respect to the resulting inner product ( , )g, Bozejko and Speicher thus 
constructed and introduced the so-called q-Fock space ^q. 

Following much interest over the course of the two subsequent decades, the g-Fock space is now 
known to have a number of remarkable properties. Far from attempting to overview them all, we 
now focus on the probabilistic aspects of the algebras of bounded linear operators and point out a 
few key results. Namely, starting with the very motivation, it was shown in [2j that the deformed 
field operator Sq{h) := aq[h)* + aq{h) is a natural deformation, in the setting of (g-CR), of the 
Gaussian random variable. The corresponding analogues of the complex Gaussian random variable 
were also constructed [TT] . More generally, an analogous formulation was found to give rise to a 
g-deformation of the Brownian motion [21 [12] and, even more broadly, to a characterization of the 
g-Levy processes [13]. Yet, there is also much about the structure of the operator algebras on the 
g-Fock spaces that remains mysterious and is object of current research. For example, remaining 
in the probabilistic vein, it was very recently shown that the g-Gaussian random variables are 
(somewhat unexpectedly) freely infinitely divisible |14j . 

The goal of this paper is to introduce a second-parameter refinement of the g-Fock space, for- 
mulated as the {q,t)-Fock space ^q,t- The (g,t)-Fock Space is constructed via a direct refinement 
of Bozejko and Speicher's framework [2], yielding the g-Fock space when t = 1. Before overviewing 

^In the sense that for any continuous function / supported on the spectrum of Si, J fdfi = ip{f{si)) with d^{x) = 
(27r)^^^^e^^ ^^dx. Equivalently, fis^) equals the n**^ moment of the standardized Gaussian random variable. 
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(g-CR) 
(1) 



the details of the construction and the main results, we take a moment to point out why the present 
refinement is, in fact, particularly natural. 



• Starting with the structural aspect and considering the combinatorial framework underly- 
ing the (7-Fock space, the notions of permutation inversions and crossings in pair partitions, 
which play the key role in the original formulation, are now replaced by the joint statis- 
tics given by permutation inversions and co-inversions and crossings and nestings in pair 
partitions (cf. Section [2] and Section |4]). In particular, while the permutation co-inversions 
are the complements of the permutation inversions, a satisfying characterization of the joint 
distribution of crossings and nestings of pair partitions is an open problem in combinatorics, 
of relevance to broader combinatorial questions \15\ [TBI, E] • 

• Whereas the distribution of the g-Gaussian operator is the unique measure that orthogonal- 
izes the g-Hermite orthogonal polynomials, given by the three-term recurrence zHn{z; q, t) = 
Hn+i{z;q,t) + [n\qHn-i{z; q,t), the {q,t)- Gaussian orthogonalizes the {q,t)-Hermite or- 
thogonal polynomials. The latter are given by the recurrence zHn{z; q, t) = Hn+i{z; q, t) + 
[n]q^tHn-i{z; q,t). (cf. Definition/Theorem |3] below.) 

• The q = < t case corresponds to a new single-parameter deformation of the full Boltzmann 
Fock space of free probability |18[ [T9] and of the corresponding semicircular operator. The 
corresponding measure is encoded, in various forms, via the Rogers- Ramanuj an continued 
fraction (e.g. j2D]), the Rogers-Ramanujan identities (e.g. [20]), the i-Airy function {21], the 
t-Catalan numbers of Carlitz-Riordan |22l [23] , and the first-order statistics of the reduced 
Wigner processes [Ml |25] . 

At this point, it is also important to note that the {q, t)-deformed framework independently arises 
in a more general asymptotic setting, through non-commutative central limit theorems and random 
matrix models. This is the subject of the companion paper |26j . In particular, analogously to the 
g-annihilation and creation operators which arise as weak limits in the Non-commutative Central 
Limit Theorem introduced by Speicher |6], the (g, t)-annihilation and creation operators appear in 
the general form of the theorem developed in [26]. In broad terms, this second-parameter refinement 
is a consequence of the passage from a commutation structure built around commutation signs, 
taking values in { — 1, 1}, to a more general structure based on commutation coefficients taking 
values in M. 

However, at this point it should also be remarked that, despite many analogies with the original 
formulation, the {q, t)-Fock space has a surprising property that distinguishes it from the g-Fock 
space when t ^ 1. Namely, the vacuum expectation state (p (cf. Section [s]) is not tracial on the 
*-algebra generated by the field operators {sg,t(^)}he^- 

The following is an overview of the main results in this paper, encompassing an overview of the 
{q, t)-Fock space construction. 

Definition/Theorem 1. Let n i— )• ui^^ denote the unitary representation of the symmetric group 

Sn on J^®" given by C/i"^ hi ■ ■ ■ hn = /i7r(i) <8> • • • "S) ^7r(n) a'^'^; for every permutation tt G Sn, 
let inv{Tr) and cinv{7r) respectively denote the inversions and co-inversions of a permutation (cf. 
Section^. Given the vector space ^ given in (Til), consider the "projection" operator Pq^t : ^ ^ ^ 





(3) 



Gonsider, further, the sesquilinear form { , )q^t on ^ given, via the usual inner product { , )q on 
the full Fock space (cf. Section^, by {rj^Clq^t = {^,Pq,t'n)o Vry,^ G 
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Then, for all n S N, P^^ is strictly positive definite for all \q\ < t, { , )q^t is an inner product, 
and the {q,t)-Fock space is the completion of ^ with respect to the norm induced by { , )q^f 

Definition/Theorem 2. Given the (q, t)-Fock space ^q^t o-nd the underlying real Hilbert space J^, 
define the respective creation and annihilation operators Uq^tif)* and ag^tif) on ^q^t, for f £ H 
and \q\ < t < 1, by linear extension of: 

aq,t{f)*^ = f, aq,t{f)*hi® ...®hn = f ®hi® ...®hn (4) 

and 

n 

Oq^tif)^ = 0, aq,t{f)hl ® ...®hn = Y^ /ife)^ hi ® . . . ® Ilk ® ■ ■ ■ ® hn- (5) 

fc=l 

Then, Oq^tif) o-nd Uq^tif)* ai"e adjoints with respect to { , )q^f Further, for all f,g £ J^, the 
operators satisfy the {q,t) -commutation relation 

aq,t{f )aq,t{gT - q aq^t{g)* aq,t{f) = (/, g) ,^ t^ , {{q, t)-CR) 

where t^ is the operator on ^q^t defined by the linear extension oft^il = Q. and t^hi . . .0 hn = 
t"'hi (g) . . . (g) /i„ for all hi, . . . , hn E J^. Moreover, for all n £ N and e(l), . . . , e(2n) G {1, *}, the 
corresponding mixed moment of the creation and annihilation operators on ^q^t is given by 

'fq,t{aq,t{hl)<^^ . . . a,,t(/l2n--l)^('"-')) = (6) 



...^{aq,t{h^^)<^-)aq,t{h,,y^^-)^, (7) 

where P2{2n) is the collection of pair partitions of[2n] (cf. Section^, eachY G P2{2n) is (uniquely) 
written as a collection of pairs {(it^i, zi), . . . , z„)} with wi < ... < Wn and Wi < Zi, with 
crossi^y) and nest{'f) denoting the numbers of crossings and nestings, respectively, inY( cf. Sec- 
tion^. Finally, for any f G , the operator aq,t{f) on JFq^t is bounded for < |g| < t < 1, with 
norm given by 



\aqAm = \\a,Afr 



.if 



tn*-qn, 

t~q 



0<-q<t<l, 
0<q<t=l 

.iT < q < t < 1 



(8) 



for 



• log(l-g)-log(l-t) 
log (t) - log {q) 



(9) 



Definition/Theorem 3. For h G the (g, t)-Gaussian element Sq^tih) G SS{^qA is given by 
Sq,t{h) := ag.t(/i) + Oq^tih)* . The {q,t)- Gaussian is self-adjoint, with moments given by 

^qASq,t{hf''-') 







<Pq,t{Sq,t{h) 



2n\ 



1 



mq,t^ 



(10) 

(11) 



1 



where [z'^]{-) denotes the coefficient of the term in the power series expansion of {■) and for all 
n E N, 



i=i ^ ^ 

Furthermore, the distribution of Sq^t{^), with \\e\\^ = 1, is the unique real probability measure that 
orthogonalizes the {q, t)-Hermite orthogonal polynomial sequence given by the recurrence 

zHniz;q,t) = Hn+i{z;q,t) + [n]g^tHn-i{z; q,t), (13) 

with Ho{z; q,t) = 1, Hi{z; q,t) = z. 

Definition/Theorem 4. For q = < t < 1, =^(=^o,t) is given as the von Neumann algebra 
generated by {aq^tih)}/^^^^'. The t-semicircular element is the q = < t < 1 specialization of the 
{q,t)- Gaussian element So,t(/i) G =^(^o,t)- The moments of the t-semicircular element are given by 

r<^NC2{2n) 

where NC2{2n) denotes the lattice of non-crossing pair-partitions and Cn ^ are referred to as the 
Carlitz-Riordan t-Catalan numbers [22l[23], given by the recurrence 

n 

C^' = Y.^'^'ctldt„ (14) 

k=l 

with Cq^ = 1. The moments of the normalized t-semicircular element so,t := so,f(e), for \\e\\^^ = 1, 
are encoded by the Rogers-Ramanujan continued fraction as 

E ^q,t{s,Ahr)z- = ^^0^ • (15) 

n>0 1 



t^Z 



1 



t^Z 



Furthermore, the Cauchy transform of the corresponding normalized t-semicircular measure /io,t is 
given by 



r^,t\ n „ z AAl/izt))' 



(l_t)(l_t2)...(l_tn) 

where At denotes the t-Airy function of [21], given by 

n2 



Mz) = yi r^-^ A-zT- (17) 

^ ^ ^ (l-t)...(l-t")^ ^ ^ ^ 



n>0 



In particular, letting {zj}j(zfq denote the sequence of zeros of the rescaled t-Airy function At{z/t), 
the measure fio^t is a discrete probability measure with atoms at 



±^t/z,, j e N 

with corresponding mass 

Mzj) 



2zjA[izj/t)' 

where A't{z) := j^At{z). The only accumulation point of no,t is the origin. 



The normalized t-semicircular measure //o,t is also the unique probability measure orthogonalizing 
the t-Chebyshev II orthogonal polynomial sequence {[/^(-z; i)}n>o> « specialization of the orthogonal 
polynomials of Al-Salam and Ismail [27] and determined by the three-term recurrence 

zUniz;t) = Un+iiz;t)+t''-^Un-i{z]t), 

with 

Uo{z;t) = 1, Ui{z;t) = z. 

Finally, t-semicircular measure /Uo,t is, in a certain sense, the weak limit of the first-order sta- 
tistics of the reduced Wigner process [HI [25]. Specifically, for all p G [0, 1] and n £ N, 

where (fjq = -^Tr(g)E and {WN,p{k)}k&i is a sequence of Wigner matrices with correlations 

EK,(A;)w,,'(^)) - I 0, otherwise ^^^^ 

¥.{wij{k)w^^j{m)) = p™-^ form>k, (20) 

with Wij{k) denoting the {i,jy^ entry of WN,p{k) . 

Finally, recall that the crux of the companion paper [26] is the extension of Speicher's Non- 
commutative Central Limit Theorem [6], giving rise to an an asymptotic model for operators sat- 
isfying the commutation relation {q,t)-CR. It should be emphasized that this asymptotic model 
provides an existence proof, independent of the explicit construction of Section [Sj and provides an 
alternative reason for the fundamental ordering bound present throughout, namely \q\ < t. 

1.1. Deformed Quantum Harmonic Oscillators. In physics, the oscillator algebra of the quan- 
tum harmonic oscillator is generated by elements {1, a, a* ,N} satisfying the canonical commutation 
relations 

[a,a*] = l, [N,a] = -a, [N,a*]=a*, 

where a* , a, and N can be identified with the creation, annihilation, and number operators on the 
Bosonic Fock space. Physicists may also speak of generalized deformed oscillator algebras, which 
are instead generated by elements satisfying the deformed relations 

a*a = f{N), aa* = f{N + 1), [N,a] = -a, [N,a*] = a*, 

where / is typically referred to as the structure function of the deformation. While an in-depth 
review of single-parameter deformations of the quantum oscillator algebra is available in [28j, of 
particular interest are the so-called Arik-Coon g-deformed oscillator algebra f29] given by 

1 - o" 

aa* -qa*a = l, [N,a] = -a,[N,a*] = a* , f(n) = —, g G M+ 

l-q 

and the Biedengarn-Macfarlane q-deformed oscillator algebra [301 [3T] given by 

aa* -qa*a = q'^, [N,a] = -a, [N,a*]=a*, /(n) = ^3^^^, ^ G M+. 

For q G [—1,1], the Hilbert space representation of the Arik-Coon algebra, generalized to an 
infinite-dimensional setting, is given by the (;-Fock space of Bozejko and Speicher [2j . Manipulating 
Definition/Theorem 1 and 2, the reader may readily verify that for q £ (0, 1], a Hilbert space real- 
ization of the Biedengarn-Macfarlane algebra is similarly given by the {q, t)-Fock space specialized 
to t = q~^. More generally, the Chakrabarti-Jagganathan oscillator algebra [32] is given by 

aa* — qa*a = p~^ , [N,a] = —a, [N,a*]=a*, f{n) = — ^ , g,pGM+. 
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For < |g| < p^^, the reader may also verify that the {q,t)-Fock space again provides the desired 
Hilbert space reahzation for for t = p~^. A Bargmann-Fock representation of this algebra was 
previously considered in [33] for the case of a single oscillator, but no explicit underlying Hilbert 
space was constructed nor shown to exist in the parameter range considered. Instead, an explicit 
representation was provided in the space of analytic functions via suitably deformed differential op- 
erators. Note that the two-parameter deformation of the Hermite orthogonal polynomial sequence 



given by (13) also appears as recurrence (15) in |33] . 

1.2. General Brownian Motion. From a high-level perspective, the ((7, t)-Gaussian processes 
fall under the framework of Generalized Brownian Motion |33|- The latter is described by families 
of self-adjoint operators G{f), where / belongs to some real Hilbert space Jif, and a state (p on the 
algebra generated by the G{f), given by 

v9(G(/i)...G(/2„))= Yi ^(^) n (/-/.•)-^^ 

for some positive definite function r : ,^^2(2?^) — ?• M. Thus, the g-Brownian motion is given 
by Tqiy^) = whereas in the present context, the ((7, t)-Brownian motion corresponds 

to Tq^ti'^) = Many other known r functions exist. Most generally, a beautiful 

framework by Gu^a and Maassen |35l |36] , proceeding via the combinatorial theory of species of 
structures [37j, encompasses the familiar deformations (Bosonic, Fermionic, free, as well as [2|l3l]). 
It is foreseeable, though not presently clear, whether (or how) the current formulation can be 
encompassed within the same framework. 

Arguably the most closely-related framework to that of the {q, t)-Gaussians is provided by 
Bozejko in [38], introduced with the goal of extending the ^-commutation relations beyond the 
q G [—1,1] parameter range. In particular, Bozejko studied operators satisfying two types of de- 
formed commutation relations, namely 

A{f)A*{g) - A*{g)A{f) = q''{f,g).^ for g > 1 

and 

B{f)B*{g) + B*{g)B{f) = \qf{f,g), for q < -1. 
Comparing Bozejko's creation and annihilation operators and his g-deformed inner product with 
the present definitions, Bozejko's setting turns out to be that of the {q, i)-Fock space for g 1— t- 1 an 
t I— 7- g. The difficulty encountered by Bozejko, pertaining to the fact that the Gaussian element 
is no longer self-adjoint, is consistent with the fact that for t > 1, the creation and annihilation 
operators are no longer bounded. 

As later found, a relevant two-parameter deformation of the canonical commutation and anti- 
commutation relations was previously studied by Bozejko and Yoshida [3^, as part of a more 
general framework. In |39], a general version of the n-dimensional "projection" operator P^^ was 
given by HILi Y^neSn Q^^^^^\ where is any sequence of positive numbers. The reader may 
readily verify that the [q, t)-projection operator P^"^ is recovered for r„ = and by substituting 
q I—)- q/t. Instead, Bozejko and Yoshida focus on the specialization r„ = s^" and q ^ q, yielding 
the so-called {q, s)-Fock space. The corresponding combinatorial structure (now given in terms of 
crossings and inner points) , the Wick-type formulas, and continued fractions are considered in [39]. 

In the context of the generalized Brownian motion, the point of the present pair of articles, which 
construct and describe various aspects of the {q, i)-Fock space, is to argue that the framework at 
hand is ultimately a highly natural refinement of the g-Fock space. This argument is based on 
the structural depth of the {q, t)-deformed framework, as evidenced by its intimate ties to various 
fundamental mathematical objects, as well as to the correspondence with the natural generalizations 
of non-commutative asymptotic frameworks (developped in [26j). 
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Figure 1. An example of a crossing [left] and nesting [right] of a pair partition 
y = {{ei,zi), {en,Zn)} of [2n]. 

o^^^^^^b^ ob^^'^^^^b^ o o'6^h^ o 
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Figure 2. Example of three pair paritions on [2n] = {1,...,6}: cross ('^) = 
3, nest(ri) = [left], cross(r2) = 2, nest(r2) = 1 [middle], cross(r3) = 
0, nesti^s) = 3 [right] 



2. Combinatorial Preliminaries 



The present section briefly overviews the key combinatorial constructs that will, in the subsequent 
sections, be liberally aplied to encode the structure of {q, t)-Fock spaces, the creation, anihilation, 
and field operators on these spaces, as well as the mechanics of the relevant non-commutative limit 
theorems. In broad terms, the objects of interest are set partitions, permutations, lattice paths, and 
certain combinatorial statistics thereof. 

2.1. Partitions. Denote by ^{n) the set of partitions of [n] := {1, . . . ,n} and by ^2(2?^) the set 
of pair partitions of [2n], that is, the set of paritions of [2n] with each part containing exactly two 
elements. (Note that the pair partitions are also referred to as pairings or perfect matchings.) It 
will be convenient to represent a pair partition as a list of ordered pairs, that is, ^2(2?^) 3 Y = 
{{wi, zi), . . . , {wn, Zn)}, where Wi < Zi for i G [n] and wi < . . . < Wn- 
Of particular interest are the following two statistics on ^2(2n). 

Definition 1 (Crossings and Nestings). For y = {{wi, zi), . . . , {wn, Zn)} G ^2(2^), pairs {wi, Zi) 
and {wj,Zj) are said to cross if Wi < Wj < Zi < Zj. The corresponding crossing is encoded by 
{wi, Wj, Zi, Zj) with 

Cross{y) := {{wi, Wj, Zi, Zj) \ {wi, Zi), {wj, Zj) G Y with Wi < Wj < zi < Zj}, 
cross{y) := \Cross{y)\. 

For y = {{wi, zi), . . . , {wn, Zn)} G ,f3^2(2n), pairs {wi, Zi) and {wj, Zj) are said to nest ifwi < Wj < 
Zj < Zi. The corresponding nesting is encoded by {wi,Wj,Zj,Zi) with 

Nest{y) := {{wi,Wj, Zj, Zi) \ {wi,Zi), {wj,Zj) G y with Wi < Wj < Zj < Zj}, 
nest{y) := \Nest{y)\. 

The two concepts are illustrated in Figures^ and^ by visualizing the pair partitions as collections 
of disjoint chords with end-points labeled (increasing from left to right) by elements in [2n]. 

Let [n]g^t denote the ((7, t)-analogue of a positive integer n, given by 

[n],,t := t"-^ + gt"-' + . . . + q^-' = (t" - <?")/(* - q). (21) 

Note that letting t = 1 yields the usual (/-analogue of integers, i.e. [n]q = riq^i. In this notation, 
both the generating functions of crossings in =5^2 (2?^) and the joint generating function of crossings 



and nestings in ^2(2n.) admit elegant continued fractions, given by 



E cross{'^) n \ ^ cross('/^) jMest(y) n 



1 _ 1 _ ["^^g'*^ 

The above continued fractions can be obtained via a classic encoding of weighted Dyck paths (see 
|40j . also in a more relevant context [H] and |15j). 

The generating function of crossings in ^2(2^-) admits an interesting explicit (even if not closed- 
form) expression, namely 

^ ^ (l - q)n 2^ ^ \n + kj' 

r&3'2{2n) ^ ^' k=-n ^ ^ 

known as the Touchard-Riordan formula |42^ H3l H5] . No analogue of the above expression is 
presently known for the joint generating function of crossings and nestings. 

2.2. Permutations. Let Sn denote the permutation group on n letters. . 

Definition 2 (Inversions and Coinversions). Given a £ Sn, for n > 2, the pair i,j £ [n] with 
i < j is an inversion in a if (j{i) > cr{j). The corresponding inversion is encoded by {i,j) and the 
set of inversions of a is given by Inv{a) := {(i,j) | i-,j S < j,a{i) > cr{j)} with cardinality 

inv{a) := \Inv{a)\. Analogously, the pair i,j G [n] with i < j is a coinversion in a if c^i) < cr{j). 
The corresponding coinversion is encoded by {i,j) and contained in the set Cinv{a) := {{i,j) \ i,j £ 
[n],i < j,cr(i) < cr(j)} with cardinality cinv{a) := \ Cinv{(T)\. For n = 1, the sets of inversions and 
coinversions are taken to be empty. 

It is well known that the generating function of the permutation inversions is given by the 
so-called g-factorial, namely 

n 

The above expression is in fact readily obtained as a product of the generating functions of the 
crossings incurred by the "chord" i 1— )• a{i) from the chords j 1— )• a{j) for j > i. By adapting this 
reasoning to coinversions, or by realizing that cinv(o") = (2) ~ iiiv(c7), the reader may verify that 

n 

r^'^'^h^'^^^-^ = i[\i]g,f (23) 

For n > 2 and o" G 5"^, it is convenient to visually represent a in a two-line notation. Then, 
a corresponds to a bipartite perfect matching and the inversions correspond to crossings in the 
diagram and the coinversions to "non-crossings". For instance, in adopting this representation, it 
becomes clear that a permutation and its inverse have the same number of inversions (and therefore 
also of coinversions). 

Remark 1. Representing a permutation cr £ Sn in a two-line notation, aligning the two rows, 
and relabeling yields a unique pair-partition vr G ^2(2^). By additionally reversing the order of 
the bottom line, as illustrated in Figure the inversions will correspond to crossings of the pair 
partition and the co-inversions to the nestings. Indeed, given an inversion i < j , a{i) > cr{j), this 
transformation yields the pairs (i, 2n-\-l — a{i)) and (j, 2n-\-l — a{j)). Since i < j < 2n + l — cy{i) < 
2n -|- 1 — the two pairs now form a crossing in n. Similarly, given a co-inversion k < m, 
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Figure 3. Example of transforming the permutation a = (231) into a pair- 
partition, with permutation inversions corresponding to crossings and coinversions 
to nestings. 

cr{k) < (T{m), the transformation yields the pairs {k,2n+l—a{k)) and {m,2n + l — a{m)), resulting 
in the ordering k < m < 2n + 1 — a{m) < 2n + 1 — cr{k) and yielding an inversion in vr. 

Naturally, this transformation of a permutation into a pair partition is by no means surjective, 
and is illustrative of the reasons why the set partitions (with the crossings and nestings) and the per- 
mutations (with inversions and coinversions) will both be found to feature in the algebraic structure 
of the {q,t)-Fock space. 

2.3. Paths. In the present context, paths will refer to finite sequences of coordinates in the lattice 
Z X Z. A North-East/South-East (NE/SE) path will refer to a path {xo,yo), {xi,yi), . . . , {xn,yn) 
where {xo,yo) = (0,0), Xi = Xi-i + 1 and y^+i e {y^+i + - 1} for i £ [n]. Indeed, inter- 

preting the coordinates as vertices and introducing an edge between (xi_i, 2/j_i) and {xi,yi) for 
i S [n] yields a trajectory of a walker in the plane, starting out at the origin and moving, at each 
step of length ^/2, either in the NE or SE direction. A Dyck path of length 2n is NE/SE path 
{xo, yo), {xi,X2), . • . , {xn, Vn) where yi > for all i G [n — 1] and yn = 0. Given n G N, the set of 
Dyck paths of length 2n will be denoted by D„ . The reader is referred to Figure [5] of Section |4] for 
an illustration. 

Counted by the Catalan number C„ = (n + 1)^^(^^), Dyck paths are found to be in bijective 
correspondence with a surprising number of combinatorial objects. In the present context, in 
Section |4j weighted Dyck paths will be found to encode the moments of the creation and field 
operators on the (g, t)-Fock space. 

3. The (g, t)-FocK Space 

The present section constructs the {q, t)-Fock space as a refinement of the construction introduced 
in [2]. 

Consider a real, separable Hilbert space Jif and some distinguished vector disjoint from J^. 
Let ^ = {CQ) © ©„>i J^^"', where J^c is the complexification of and both the direct sum and 
tensor product are understood to be algebraic. In particular, ^ can be viewed as the vector space 
over C generated by {Q} U {/ii . . . /in}/iie.^,neN- 

For / € and q, t G M, define the operators a{f)* and a{f) on ^ by linear extension of: 

a(/)*0 = /, a{f)*hi® ■ . .(^hn = f ®hi(^ . . .®hn (24) 

and 

n 

a{f)n = 0, a{f)hi(g) ...(g)hn = ^q'''^t"^'''{f,hk)^hi(g) ...(g)hk(g) ...(g)hn, (25) 

k=l 

where the superscript hk indicates that h^ has been deleted from the product. Note that the 
creation operator a(/)* is defined identically to the operator c*(/) in [2], whereas the "twisted" 
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annihilation operator a(f) is the refinement of the operator c(/) in [2] by a second parameter, t. 
Next, define the operator on ^ by hnear extension of 

t^n = n, t^hi ® ...®hn = f'hi® ...®hn. (26) 

For t > 0, the operator can be written in a somewhat more natural form as e"^, where a = log(t). 

Lemma 1. For all f,g €z .3^ , the operators a{f),a{g)* on ^ fulfill the relation 

aUM9r-qa{graif) = {f,g)^t^. 



Proof. For any n G N and g,hi, . . . ,hn G J^, 

a{f)a{g)*hi ^ . . . hn = a{f)g hi (g) . . . hn 

n+l 

= g)^ hi(g)...(g)hn + Y^ hk^i)^g ®hi® h-i O . . . O /in 

k=2 

= g)^ /ii «)...«) /i„ + 5 «) q^f'-^if, hk),^ hi0 ...0hk^ ...0hr^ 



f'if, g)^hi(^...(^hn + q a{g)*a{f)hi Cg . . . hn 



□ 



Define the sesquilinear form (,)<j,t on ^ by 

{gi ■ . . (S) gn,hi(g) . . . hm)q,t = for m / n (27) 
and otherwise recursively by 

n 

{gi®...®gn,hi®...®hn)q,t = ^q^'^ ^'-^'{gi.hk) ^{g2®- ■ -^gn^hi®. . .®hk®. . .®hn)q,t- (28) 

fc=i 

Remark 2. For all q, setting t = 1 recovers the usual inner product on the q-Fock space of [2]. 
Letting t = yields a sesquilinear form on the full Fock space that is given by 

{gi . . . gn,hi ^ . . . hn)qfl = q^^\gi,hn)^ ■ • • {gn,hl)^, 

which does not generally satisfy the positivity requirement of an inner product. 

The range of t required for (, )q^t to be an inner product will be characterized shortly, in Lemma|4j 
But, first, note that for all / € J^, a{f)* is indeed the adjoint of a{f) with respect to (, )g^f 

Lemma 2. For all f e M", ^,ri e ^, 

{a{frC,v)q,t = {ta{f)v)q,t. 

Proof. It suffices to note that, directly from the previous definitions, 

{a{f)*gi <^ . . . <^ gn,hi0 . . . hn+i)q,t = {f (g gi ® . ■ ■ (g gn,hi <S) . . . ® hn+l)q,t 

n+l 

^ ^ ^ ^n+l-k^j^ f^^^^ {gi0 ...0gn,hi(g) ...^hk(S> ...^ K+l)q,t 

k=l 

= {gi<^ ■ ■ gn, a{f)hi ® hn+l)q,t 

□ 
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Still in line with [2], it is convenient define the "projection" operator Pq^t : ^ ^ ^ allowing 
one to express the sesquilinear form (, via the usual scalar product (, )o on the full Fock space. 

Consider the unitary representation vr i— )• C/i-"'' of the symmetric group Sn on given by 

Recalling the permutation statistics given by inversions and co-inversions, defined in the previous 
section (cf. Definition [2]) , let 

Pl,t = ^Pif with Pg(? : ^ 

ra=0 

:= qinvin)^cmv(n)jj{^n)^ (29) 

where for the unique vr E 5*1, it is understood that inv(7r) = cinv(7r) = 0. Thus, Pg^J^ = 1 and 
there is no change to the scalar product on the single-particle space. For n > 2, note that for 
every 1 < i < j < n, the pair (i, j,7r{i),TT(j)) is either an inversion or a coinversion, and therefore 
inv(7r)+cinv(7r) = (2). It follows that for n >2 and t 7^ 0, the projectiorj^operator can equivalently 
be expressed as 

pw = ,ra 2: (?)'■"'"' f/<»> = <R)F<;>. (30) 

where Pg"^ denotes the projection operator on the subspace of the g-Fock space (see [2j). 

Lemma 3. For all ^, G 

{V,Oq,t = {tPq,tV)0- 



Proof. It suffices to prove that for all n G N and all gi, hj £ J^, 

{gi® ...®gn,hi® ...®hn)q,t = {gi® gn,Pq,thl® ...®hn)Q 

The claim clearly holds for n = 1. Proceding by induction on n, recall that from the definition of 

(> )q,t, 

n 

(51 . . . O . . . hn)q,t = ^g''"-^r~''(5'l,/ifc)^(52 (g) . . . (8)5n,/il ® • • • (8) /ifc . . . /ln)g,t- 

fc=l 

Letting s!^^ denote all bijections from {2, . . . , n} to {1, . . . ,k, . . . , n}, note that, for all n E N, every 
n £ Sn can be uniquely decomposed as a pair (A;, a) for some k G [n] and a G 'S'^^;^ and, conversely, 
that any such pair gives a distinct element of 5„. Specifically, let 7r(l) = k and 7r(£) = a{i) 
for i G {2, . . . ,n}. Then, noting that the natural correspondences [n — 1] -H- {2, and 
[n — 1] -H- {1, . . . , A;, . . . n} are order-preserving, the inductive hypothesis on n — 1 can be written as 

{g2 ...0gn,hi(^ ...0hk(^ ---(^ hn)g,t = Yl <i'"^''^t^''"^''H92, K{2)).^ ■ ■ ■ {9n, K(n)).^, 



■^Strictly speaking, Pq,t is not a projection for general values of q and t, as ^-P,"'^ = X^irss (/■^'"^''^^f^'^'"^''^' Ui"'\ 
which follows from the fact that a permutation and its inverse share the same number of (co)inversions (cf. Section[2|. 
The terminology is inherited from the classical construction, as Pq^t = Pq,t if (and only if) |g| = f = 1. 
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1 n 

QO O - O- O O 

6 6 6 - o- o o 

V . ^ 

k - \ additional n - k additional 
inversions coinversions 

Figure 4. If 7r(l) = fc, then there are exactly k — \ elements from the set {2, . . . , n} that 
map under tt to an element in {1, . . . , — 1} and exactly n — k elements from {2, . . . ,n} 
that map to {fc + 1, . . . , n}. Thus, there are k — 1 elements j > 1 for which < 7r(l) and 
n — k elements ^ > 1 for which Tr{£) > tt{1). For any other pair € [n]^ with i ^ 1, the 
corresponding inversion or coinversion is given by the map a : {2, . . . , rt} — > {1, . . . ,k, . . . , n}. 
It follows that inv(7r) = fc — 1 + inv(f7) and cinv(7r) ^ n — k + cinv((T). 



where inv(o") counts all the pairs (i, j) S {2, . . . , n} x {1, . . . , k, . . . , n} with i < j and 7r(f) > 7r(j) 
and cinv(a") is defined analogously. Furthermore, observing that inv(7r) = inv(cj) + A; — 1 and 
cinv(7r) = cinv((T) + n — k, as demonstrated in the caption of Figure |4| it follows that 

(51 (g) . . . (g) 5„, /ll (g) . . . (g) hn)g^t 
n 

□ 

While the following facts, contained in Lemmas |4] through [5} have direct proofs analogous to 
those in [2j , it is more convenient to use ( 30 ) and derive the desired properties from those of the 
g-Fock space. 

Lemma 4. a) The operator Pq^t is positive for all \q\ < t. 
b) The operator Pq^t is strictly positive for all \q\ < t. 

Proof. Since Pq^t = ©.^o-^g"^' suffices to consider the positivity of Pq^^ ■ Since t^a) > 0, the 
positivity of -P^"^ for |g| < t follows from (30) by the positivity ([2j) of pj^J^ for \q/t\ < 1 and the 
strict positivity of of P^"^ for Ig] < t follows from the strict positivity ([2J) of P^jj for \q/t\ < 1. □ 

Since ( , )q^t is an inner product on the completion of ^ yields the desired {q, t)-Fock space. 

Definition 3. The {q,t)-Fock space ^q^t is the completion of ^ with respect to { , )q^f 

Remark 3. For dim{M') > 2, the conditions of Lemma are also necessary. In particular, let 
61,62 be two unit vectors in with (61,62)5,4 = 0. Note that 

||6i (g 62 + 62 (g eillg^f = 2t + 2g, and ||6i (g 62 - 62 (g 6i||q,t = 2t - 2g. 

Thus, { , )q^t is positive (resp. strictly positive) only if \q\ < t (resp. \q\ < t). 

Analogously the case of the g-Fock space of [2j, the operators a{f) on ^g^t are bounded for 
< —q < t < 1 and < q < t < 1. Letting q = 1 (and therefore t = 1) recovers the Bosonic Fock 
space, in which case a(/)* and a{f) are unbounded and defined only on the dense domain 

13 



Lemma 5. For any f S the operator a(f) on ^ is bounded for < |g| < t < 1, with norm 
given by 

U 0<-q<t<l, 



Hf)\\ = 



n-q 



fn* -qn* 

t-q 



< < t = 1 



< (7 < t < 1 



for 



• log(l-g)-log(l-t) 
log (t) - log (g) 



Proof. Let q £ [-1, 0] and \q\ < t < 1. For € Jf®", by Lemma[T| 

{aifYCn, aifrCn)q,t = {in, a{f)aifrCn)q,t = (/, f)j^ (Cn, ^n),,* + qi^n, aiff a[f)in) q,t 
= 11/11^ llen||^, + g||a(/Kn|lL<||/|l3r'llen||^„ 

as g < and < t < 1. Next, for an element ^ = Y17=o ^ where G C and .^j G 
that a{f)* is linear and {(,m£,m)q,t = whenever n ^ m implies that 

n n 

{a{frta{froq,t = Ei"^i'w)*e^,a(/n.)a,* < = ii/ii^iieiiL- 

i=0 i=0 

Finally, since (a(/)*f], a(/)*17)g,t = \\f\\%, it follows that ||a(/)|| = ||a(/)*|| = ||/||^. 

For q £ (0,1), analogously to the previous case, it suffices to focus on ^ G J^®". Again, 
{aifrn,a{frn)g,t = \\f\\%' and, for n G N, Jsof yields 

Hfn, aifroq,t = {aifTt p^T'^aifTOo = t("t') (a(/)*e, p^;+'^a(/)*e>o 

fn + l\ 



t^"^^'Haifri,aifrOq/t- 



Recalling that the creation operators on the (g, t)-Fock and /x-Fock spaces (for fi G [—1,1]) are 
defined by the linear extension of the same operator on the dense linear subspace it follows that 

{a{fr(,a{frOq/t = {aq/t{fri,a,/t{frOq/t, 
where a*^^{f) analogously denotes the creation operator on the (g/t)-Fock space. Recalling, fur- 
thermore, the familiar bounds on ||aq(/)Cllg (s-g- Lemma 4 in |2|), 



But, 



and so. 



2 



{a{frLaifrOq,t<t' 



n+1 



Mf\\%m\ 



Since < t < 1, it follows that a(/)* is bounded. To recover the corresponding expression for the 
norm, let ^ = / O . . . / = G ^'^^ and therefore a{f)*C = Thus, 



%f\\i\\l,f 



n+l 



{a{fri,a{froq,t = = j^g^-^^+i 

fc=i 
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j.n+1 _ „n+l 

-^"^".^iir"iiL= 11/113.-11^11 

I q 



q,t 



For t = 1, 



sup 

neN 



(a(/)*/®",a(/)V®")g,, 



1 



2 



Otherwise, it remains to compute n G N that maximizes t"^^ — q"^^. Let 

^n— 1 

rn-= , nGNUiO}, 

and note that, since < (7 < t, both the numerator and the denominator are strictly positive. 
Furthermore, {r„}„,>o forms a strictly increasing sequence as 



(1 - t) 



t 



-1 > Tn-l- 



(1 - g) q 

Now, note that r„ > 1 iff — (7" < t"^^ — g"^^ and, conversely, r„ < 1 iff — > — g""^. 
It follows that t" — is maximized for n = n^,, where n^, is the greatest non- negative integer for 
which r„ < 1. A straightforward calculation then yields 

"log(l-g)-log(l-i)" 



log (t) - log {q) 



□ 



4. The (g', t)-GAussiAN Processes 

The natural starting point to the probabilistic considerations of this section is the *-probability 
space {^q^t,^q,t)i discussed next, formed by the *-algebra generated by the creation and anni- 
hilation operators on ,!^q^t and the vacuum-state expectation on the algebra. The *-probability 
space {'^q,tiVq,t) provides a convenient setting in which the (g, t)-Gaussian family is subsequently 
introduced and studied. 

4.1. The {^q^t^^q,t) *-probability space. For a clear introduction to non-commutative proba- 
bility spaces, the reader is referred to the monograph [46j. In the present context, the setting of 
interest is that of the (^^q,t, ^q,t) *-probability space, formed by: 

• the unital *-algebra ^^^.^ generated by {a{h) \ h G 

• the unital linear functional : 'iSq^t — >■ C, 6 1— {^},bil.)q^t (i-e. the vacuum expectation 
state on 5^g,t.). 

In this non-commutative setting, random variables are understood to be the elements of ^q^t 
and of particular interest are their joint mixed moments, i.e. expressions . . . 5^^^^) for 

61, . . . , 6fc G "^^q^t, k £ N and e(l), . . . , e{k) G {1, *}. Since tpq^t is linear, one is allowed to focus on 
the joint moments fq^tidihiY^^^ . . . aQikY^^'^). These turn out to be intimately connected with two 
combinatorial objects, discussed next. 

For n G N, consider first the map ipn from {1, *}" to the set of all NE/SE paths of length n (cf. 
Section [2]), by which every "*" maps to a NE step and every "1" maps to a SE step. Clearly, ipn is 
a bijection and of particular interest is the set ip2ni-^n), where denotes the set of Dyck paths 
of length 2n. An example of an element in {1, *}^^ mapping into -D7 is shown in Figure [5] 

Given (e(l), . . . , e(2n)) G let each * encode an "opening" and each 1 encode a "clo- 

sure" and, given some such string of openings and closures, consider all the ways in which the 
elements of the string can be organized into disjoint pairs so that (1) an opening is always to 
the left of the corresponding closure and (2) no opening is left unpaired. Clearly, each fixed 
(e(l), . . . , e(2n)) G il)~^{Dn) is thus associated with a distinct set of pair-partitions of [2n], that 
is, a subset of =^2(2?^) (cf. Section [2]). Conversely, every pair partition naturally corresponds 
to a unique string of openings/closures belonging to ip~^{Dn)- Thus, writing Y "V' for any 
-y ^-y £ ^2(2^) for which the corresponding strings of openings and closures both encode the 
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Figure 5. tpui*, *, *, *, *, 1, 1, *, 1, 1, 1, 1, *, 1). 



same string in ^p ^{Dn) defines an equivalence relation. Figure [g] shows the six pairings in the 
equivalence class of the string (*, *, *, 1, 1, 1) E iP~^{Dq). 



oo^^^b^^h\) o6^^^^^^ od^^^^b^h^ 



1 



1 



o o o o o o o o 

1 6 1 6 





Figure 6. The equivalence class of (*,*,*, 1,1,1) e i) (-De) in -P2(6). 



Next, recalling that the purpose is calculating the value of the mixed moment 

(/.,,i(a(/iir«...a(/i2nr(2n))^ 

fix a string (e(l), . . . , e(2n)) G ■ip2ni^n) and jointly consider some underlying choice of hi, . . . , h2n S 
J^. Each pairing y = {{w\,zi), . . . , {wn, Zn)} in the equivalence class of (e(l), . . . , e(2n)) is then 
assigned the weight 

n 

Wt(r; hi,..., h2n) ■■= l[{h^^,h,^),^, 

i=l 

where cross('^) and nest('^) correspond, respectively, to the numbers of crossings and nestings in 
Y, as defined in Section [2] For example, the top left-most pairing in Figure [6] is thus given the 
weight 9^(51, 54)^(52, fl'5>.ir(93,56).if'- For n £ N, let Tg^t{hi, . . . ,h2n;e{l), ■ ■ ■ ,€{2n)) denote the 
generating function of the weighted pair-partitions of [2n] , namely 

rg,f(/ii, . . . ,/i2n;e(l), • • • ,e(2n)) := ^ wt{y; hi, . . . ,h2n) 

r~p(£(l),...,e(2n)) 

n 

re.f?'2(2n) i=l 
r~p(£(l),...,€(2n)) 

where Y ~p (^(1)) • • • > ^C^^)) is meant to indicate that Y is in the equivalence class of (e(l), . . . , e(2n)) 
under ~p (i.e. Y has its opening/closure string given by (e(l), . . . , e(2n)) ). Writing the weight of a 
pairing as a product of weights of its pairs, as in the following lemma, is the remaining ingredient 
in connecting the combinatorial structures at hand to the moments (pq^t{0'{hiY^^^ ■ . . a(^2n)^^^"'^)- 
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1 2 3 4 A Y— ' ' V ' A 

^ — Y — ' ^ — Y — ' generate generate 

available available nestings crossings 

Figure 7. The sequential pairing of openings and closures corresponding to the 
Dyck path of Figure [5j with an arrow denoting the currently considered closure. In 
the left figure, there are four available (i.e. as of yet unpaired) openings. In the 
right figure, the current closure is paired to the indicated opening, thus incurring 
two crossings and one nesting. Note that the crossings and nestings incurred by 
the current closure never include any crossings or nestings already counted in the 
previous closure. 

Lemma 6. For Y = {{zi,wi), . . . , {zn,Wn)} G ^2(2n) with wi < W2 < ■ ■ ■ < Wn (i-e. with pairs 
indexed in the increasing order of closures) and hi, ... , /i2n G 

n 

Wt{y- hi,..., h2n) = n gl^^^l^'<^^<"'-^>*>l tl^^^|l^"^<"-^'>*>l(/l.,, V),^. 
i=l 

Proof. Consider the following procedure for assigning weights to a pairing. Starting with the left- 
most closure wi, suppose that it connects to some opening zi (where zi <wi). Since wi is indeed 
the left-most closure, exactly wi — zi pairs will cross the pair {zi,wi) and exactly zi — 1 pairs 
will nest with it. Thus, assign the pair {zi,wi) the weight q^^~^^t^^~^{hzi, hw^),^. Consider now 
the z*^ closure from the left, Wi, and let be the pairing from which all pairs {zk,Wk) for k < i 
have been removed. Then, assigning the pair {zi,Wi) its weight in 1^' according to the previous 
recipe does not take into acount any crossings or nestings that have already been accounted for 
by the previous pairs. It now suffices to note that, by the end of the procedure, all the crossings 
and nestings have been taken into account and that the product of the weights of the pairs indeed 



equals the expression for wt(';^; hi, ... , /i2n) in (31 ). □ 



(/.g,i(a(/iirW...a(/i2„)^('")) 



The sequential procedure for assigning weights to pairings is illustrated in Figure [7j We are now 
ready to calculate the joint mixed moments of interest. 

Lemma 7. For allnGN and e(l), . . . , e(2n) G {1, 
<^,,t(a(/iir«...a(/i2„-ir(=^""^)) = 

0, (6(2n),...,e(l))0V2~nPn) 

Tg,t{h2n, hi;e{2n), e(l)), (e(2n), . . . , e(l)) G i^^niDn) 
= Yl V>q,Mhn,,r'^'"'^a{h,,)<'^^) . . . <^,,t(a(/i^J<'"")a(/i,J^(^"))g-°-(^^)t"-*(^), 

where each 'Y is (uniquely) written as a collection of pairs {{wi, zi), . . . , {wn, Zn)} with wi < . . . < 
Wn and Wi < Zi. 

Proof. Given the mixed moment ipq^ticLihiY^^^ . . . a{hkY^^^), consider the reverse string e{k), e{k — 
1), . . . , e(l) and the corresponding NE/SE path ■0fc(e(/c), e{k — 1), ... , e(l)). Since a{h)Q = for all 
h £ J^, by recursively expanding the mixed moment ipg^tio-ihiY^^'^ . . . a{hkY^^^) via (24) and (25), 



it immediately follows that the moment is zero if for any 1 < i < k, the number of SE steps in 
'ijji{e{k) , e{k — 1), . . . ,e(i)) exceeds the corresponding number of NE steps. Moreover, since for all 
hi,gi G J^, {hi . . .<Sihn, gi® . . .0 gm)q,t = whenever n / m, it follows that the moment vanishes 
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unless the total number of NE steps in ipk{e{k), . . . , e(l)) equals the corresponding number of SE 
steps. This shows that the moment ipq^tiaihiY^^^ . . . a{hkY^^^) vanishes if ^/'fc(e(A;), . . . , e(l)) is not 
a Dyck path. In particular, the mixed moment vanishes if k is odd. 

Now let k be even with ij){e{k), . . . , e(l)) € Dn- For some m > 1, let e(m + 1) correspond to the 
first * from the left in the reverse string (e(A;), . . . , e(l)) (that is, corresponding to the right-most 



creation operator in a{hiY^^^ ■ ■ ■ a{hkY^^^)- Then, by (24), a{hm+i) acts on hm® ■ . .(8>/ii to produce 



a weighted sum of (m — l)-dimensional products, that is, 

m 

a{hm+i)hm ® ...®hi = q'^~''f~^{hm+i,h.i),^^h.m ® . . .hi ® ■ ■ ■ ^ hi. 

i=l 

At the same time, diagramatically, e(m + 1) corresponds to the the first closure from the left in 
ip{e{k), . . . , e(l)). Supposing that this closure pairs to the i^^ opening (from the left), for 1 < i < m, 
the weight of the resulting pair in the sense of Lemma [g] is then given by q^~'^f'~^ {hm+i,hi) 
Furthermore, the act of removing hi from the product h^ ® . ■ ■ ® hi diagramatically corresponds 
to removing the previously completed pairs in the procedure of Lemma [6] and, in both cases, the 
same iteration is subsequently repeated on the thus reduced object. Now, by definition, summing 
the weights wt(';^; hk, . . ■ ,hi) over all pairings ^ ~p (e(^)i • • • ; e(l)) is equivalent to summing the 
products of the weights of the individual pairs over all the possible ways of matching all closures 
to openings (and thus, in the above notation, over all choices of i and analogous choices made on 
the subsequent iterations). Thus, 

wt{r;hk,...,hi) 

T^p{e{k),...,e{l)) 

is exactly the sum of weights obtained by unfolding the expression (pq^tio-ihiY^^^ • • • a{hkY^^^) ^'^^ 
the recursive definitions (24) and (25). In other words, we have shown that 

^qM^iT^^^ ■ ■ ■ a(/i2n)^('"^) = Tq4h2n, • • • , hi' e{2n) , e(l)) 
whenever (e(2n), . . . , e(l)) G ip2niDn)- 

Finally, that ipg^t{a{hiY^^^ . . . a(/i2n)''''^"-*) also equals 

^ <^g,t(a(/i^J^("^^a(/i,J^(^i)) . . . (^,,i(a( V)<-")a(/i,J^(^")) g-oss(r)^nest(r) 

follows immediately from the fact that (pq^t{o.{hwjY^^^^ '^i^ZiY^^''^) — unless e{wj) = 1 and e{zj) = 
*; in other words, unless Y ~p (e(2n), . . . ,e(l)). 

□ 

In particular, the mixed moments of the single element a{h) are of a particularly insightful form. 

Lemma 8. Given a Dyck path ^„(e(2n), . . . , e(l)) and h G .3^ , let wt{€{2n), . . . , e(l); h) denote the 
weight of the path taken as the product of the weights of the individual steps, with each NE step 
assigned unit weight and each SE step falling from height m to height m — 1 assigned weight 

™ Mm m 

H.,t:=E'?""*"' = ^3^- (32) 
i=i ^ ^ 

(See the illustration of Figure^) Then, Lpq^t{a{hY^^^ ■ ■ ■ a{hY^'^'^'^) = wt{e{2n), . . . ,e{l);h). 

Proof. Returning to the sequential procedure in the proof of Lemma |6} note that when hi = . . . = 
h2n = h, the choice of an opening for any given closure only affects the weight of the present pair 
and does not affect that of the subsequently considered pairs. In particular, for any given closure 
and m available (previously unpaired) openings to the left of it, the generating function of the 
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Figure 8. Weighted path corresponding to Figure [s] for \\h\\_%- = 1. 



weight of current pair is therefore given by X^^^^ (f^ ^ll^llj^"- Moreover, by the same token, the 
sum of weights of all the pairings Y in the equivalence class of a fixed (e(2n), . . . , e(l)) G ip'^lDn) 
is given by the product of the corresponding generating functions. Finally, note that for any given 
closure, the "m" is determined by the underlying (e(2n), . . . , e(l)); specifically, the reader may 
readily verify that m is exactly the height of the Dyck path preceding the corresponding given SE 
step encoding the given closure. □ 

Figure [8] provides an example of assigning weights to the steps of the Dyck path given by 
ipi4{*, *, *, *, *, 1, 1, *, 1, 1, 1, 1, *, 1) according to the rules defined in the above proof. In particular, 
in light of Lemma [8j the product of the weights of the individual steps then yields the moment 

^,,MhHhr{a{h))Mhr{a{h))\a{hrf). 

At this point it is also interesting to note the combinatorial statistics featuring prominently in 
the previous section, given by inversions and coinversions of permutations, are intimately related to 
those of crossings and nestings in pair partitions considered presently. In particular, by Lemma [Sj 

{gnC^...C^gi,hi(E)...C^ hn)q,t = ^'''"^"^^"''"^"^^n, • • • (gl, K{n)).)^. 

while, at the same time, 

{gn® ■■■®gi,hi® ...®hn)q,t = a{gi) . . . a{gn)a* (hi) . . . a{hn)^) g^t 

n 

Y: q--'~'''h-'^'''^ll{g^^,h.^)j,, 

where the second equality follows by Lemma [Tj In other words, the value of {gn (E> ■ ■ ■ (E> gi,hi ^ 
. . . (8) hn)q^t is computed by summing the weights over pairings 'f in the equivalence class of 
(1, 1, . . . , 1, *, *, . . . , *). Writing each such Y in the two-line notation of the previous section, with 
the e(i) = 1 corresponding to the top row and e(j) = * to the bottom row, the reader may readily 
verify that the crossings turn into inversions and nestings into coinversions, and vice- versa. 

4.2. Joint and marginal statistics of {q, t)-Gaussian elements. Analogously to both classical 
and free probability and, more generally, to the non-commutative probability constructed over the 
g-Fock space of [2j, the so-called "Gaussian elements" will turn out to occupy a fundamental role 
in the non-commutative probability built over the {q, t)-Fock space. More precisely, for h G J^, the 
{q, t)-Gaussian element Sq^t{h) is defined as the (self-adjoint) element of '^^q^t given by a{h)+a{h)* . In 
[26], the (g, t)-Gaussian family will be tied to an extended non-commutative Central Limit Theorem 
and, in particular, will be identified with the limit of certain normalized sums of random matrix 
models. Meanwhile, note that the joint moments of the [q, t)-Gaussian elements come in a form 
that is the (g, t)-analogue of the familiar Wick formula (e.g. [36]). 

19 



Definition 4. For h £ .3^, let the (g, t)-Gaussian element Sq^t{h) G ^^<j,f be given by Sq^t{h) '■= 
a{h) + a{h)*. 

Lemma 9 (The q, t-Wick Formula). For a// n E N and hi, ... , /i2n G 

^q,t{Sq,t{hl) . . .Sq^t{h2n~l)) = 

^,,t{s<iAhi)---s,,t{h2n)) = q''°''^^''^t^'''^''^fl{h^^,h,^)^. 

r={{«>i,2i),...(w„,2;„)}G.5?'2(2n) i=l 



Proof. It suffices to note that 

'PgASgAhl)---Sq,t{hn))= 99,,^ (a(/li ^ « . . . a(/l2„r (^n) ) ^ 

e(l),...,e(n)G{l,*} 

from which the result follows directly by Lemma [7] and (31). □ 

Focusing on a single element Sq^tih), the corresponding moments can be expressed in a par- 
ticularly elegant form. In particular, in light of Lemma [8 the moment ipq^t{sq,t{h) ■ ■ ■ Sq^t{h)) is 
expressible via generating functions of weighted Dyck paths, which one can then interpret as a 
continued fraction via a well-known correspondence (e.g. 110]). Note that the continued- fraction 
formulation in the following lemma already features in [15j , in a more general combinatorial context. 

Lemma 10. The moments of the {q,t) -Gaussian element Sq^t{h) are given by 



2n r n] 



[1] 



where [-z"](-) denotes the coefficient of the term in the power series expansion of (•). 
Proof. 

Expanding the moment and applying Lemma [8j 

^qAsqAh?"") = Yl wt(e(2n),...,e(l);/i,...,/i) 

e(l),...,£(n)e{l,*} 

E wt{e{k),...,e{l);h). 

{e{k),...,e{l))&,lj-^Dn) 

The continued fraction expansion of J2n>o ^q,tisq,t{h)'^^)z"- is then obtained by the classic encoding 
of weighted Dyck paths (see [40j, also in a more relevant context [41] and [15J). □ 

By the previous lemma, the ffist four even moments of Sg,t(/i) are thus given by 

fq,t{Sq,t{hf) = 1 

fqA^qAhf) = \\h\\l^ 

^qASqAh)') = \\h\\%{l + q + t) 

H^qAsqAhf) = \\ht^{l + 2q + 2t + 2qt + q^ +t^ + 2qH + 2qt^ +q^ + t^), 

may be verified by enumerating all chord- crossing diagrams with, respectively, between and 3 
chords and counting the corresponding crossings and nestings. 
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Remark 4. Choosing ||e||,^ = 1, ipg^tisq,t{S')'^^) is equal to the generating function of crossings and 
nestings in i3^2(2n), that is, 

whose continued-fraction encoding, as previously noted, already explicitly features in [15] . 

Since [n]q^t = the continued-fraction expansion shows that the joint generating function is 

in fact symmetric in the two variables, that is, 

Given the fundamental constraint \q\ < t, necessary for the positivity of the sesquilinear form {,)q^t, 
this symmetry is surprising. Indeed, let 

1 n odd 

and note that associating Tq^n) with the moments of a self-adjoint element (viz. Sq^e)) in a C* 
probability space ensures the positivity of {Tq^t{n))nefi (as a number sequence) in the range \q\ < t 
|47j ■ The above symmetry then yields the positivity ofTt^q{n) for \q\ < t, which, in turn, guarantees 
the existence of a real measure whose moments are given by Tt^q{n) |47J. This observation naturally 
beckons the construction of spaces that would give rise to such "reflected" Gaussian algebras. 

Remark 5. Letting t = 1 in the moment formulation of Lemma [T^ yields an analogous contin- 
ued fraction with [n]q playing the role of the coefficients [n]q^f What is more, the corresponding 
generating function admits an interesting explicit (even if not closed-form) expression, namely 

y(^^2{2n) ^ ^' k=-n ^ 

known as the Touchard-Riordan formula 132 



The general form of continued fraction appearing in Lemma [To] is known as the Stieltjes- fraction 
(S'-fraction). Specifically, an 5- fraction is written as + + + + • • ., where, in the 
present setting, = [n]q^t. In the combinatorial approach to orthogonal polynomials pioneered 
by Flajolet and Viennot [lUl HH], an 5-fraction ^ + + + + . . . and the orthogonal 
polynomial sequence {yn{z)}n>o given by yo{z) = 1, yi{z) = 1, and 2/„+i(z) = zyn{z) - \nyn-i{z) 
can both be considered as encodings of a certain object. Namely, they both encode a collection 
of Dyck paths with NE steps carrying unit weight and each SE step of height n i— )■ n — 1 carrying 
weight \n. Returning to the setting at hand, consider therefore the following (further) deformation 
of the classical and quantum Hermite orthogonal polynomials. 

Definition 5. The {q,t)-Hermite orthogonal polynomial sequence {Hn{z;q,t)}n>o is determined 
by the following three-term recurrence: 

zHn{z]q,t) = Hn+i{z;q,t) + [n]q^tHn-i{z; q,t), 

with Ho{z; q,t) = 1, Hi{z] q,t) = z. 

While the following Lemma [IT] can also be deduced from the classical theory (e.g. [47] ) , the 
combinatorial approach of [H] reveals the rich structure associated with orthogonal polynomial 
sequences and provides a lucid link between continued fractions, moment sequences, and orthog- 
onal polynomials. In particular, interpreting the underlying Hankel determinants as certain sums 
of weighted paths, the S-type fraction + + + + . . . is readily shown to be the 
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continued-fraction encoding of the generating function X^„>o "^n-z", where {mn)n>o is the moment 
sequence associated with an orthogonahzing hnear functional for the polynomial sequence yoiz) = 1, 
yi{z) = 1, = zyn{z) — Xnyn-i{z)- In the present context, in the light of Lemma [Tof this 

correspondence yields that the {q, t)-Hermite polynomial sequence is orthogonal with respect to 
the linear functional given as the integral against the measure of Sg^4(e). In turn, since Sg^i(e) is 
bounded for \q\ < t < 1 (as the real part of a bounded operator, cf. Lemma [5]), the corresponding 
measure is compactly supported and is therefore uniquely determined by its moment sequence. 
To summarize: 

Lemma 11. The distribution of the {q,t)- Gaussian element Sq^t{e), where \\e\\j^ = 1, is the unique 
real measure that orthogonalizes the {q,t)-Hermite orthogonal polynomial sequence. 

4.3. The (Tq^t,fq,t) *-probability space. While the pair {'^^q^t,fq,t) provides a natural setting in 
which one can discuss the mixed moments of the creation operators and field operators, it is not 
an especially well-behaved non-commutative probability space. In particular, the functional ipg^t 
is neither faithful nor tracial, as (/9g^f (a(/i)*a(/i)) = ^ ipq^t{cL{h)a{h)*) = \\f\\,^. As in the full 
Fock space and, more generally, the g-Fock space, one may wish instead to focus on a subalgebra 
of non-commutative random variables generated by the (g, t)-Gaussian elements of Definition |4j 
In particular, consider the unital *-algebra generated by {sq^t{h) \ h £ J^} and let Tq^t denote its 
weak closure in I^{^q^t)- In other words, Tq^t is the ((7, t)-Gaussian von Neumann algebra. While 
the vacuum expectation state on the corresponding algebra on the g-Fock space, corresponding to 
the t = 1 case in the present setting, is tracial [3] , the same does not hold in the more general case. 

Proposition 1. For dim{Jif) > 2, Tq^t is tracial if and only ift = l. 

Proof. The forward direction is established in pj in a more general setting. For the converse, it 
suffices to consider some four vectors hi, h2, h^, E Using Lemma [9| 

^q,t{Sq,t{hl)Sq^t{h2)Sq,t{hz)Sq,t{hA))q,t = {hi, /l2 ) .if- ( /^S , /l4> J^+i (/ll , M^' (^2 , hz).^+q{hl, hz)M'{h2, hi) 

and 

^q,t{Sq,t{hi)Sq^t{hl)Sq^t{h2)Sq,t{h-i))q,t = (/14, (/i2, h^i) _^^+t{hi, h-^) ,%'{hi, h2) .^•+q{hi, h2),^{hi,h^) 

Since is a real Hilbert space, the terms in q are equal in both of the above expressions, but the 
terms in 1 and t are interchanged. Denoting by {cj} the basis of and taking hi = h2 = ei and 
/13 = /i4 = 62 , it follows that the two expressions are equal if and only if t = 1 . □ 

Remark 6. The fact that the term in q in the previous example remained the same for both 
{Sq^t{hi)sq^t{h2)sq,t{h'i)sq^t{h4)^)q,t and {VL, Sq^t{hA)Sq^t{hi)sq^t{h2)sq,t{h:i)^) q,t is not Coincidental. 
Indeed, from a combinatorial viewpoint, commuting two products of Gaussian elements is equivalent 
to rotating (by a fixed number of positions ) the chord diagrams corresponding to each of the non- 
vanishing products of the underlying creation and annihilation operators. The observed equality then 
follows from the fact that the crossings in a chord diagram are preserved under diagram rotations. 
The same is however not true of nestings, which is the combinatorial reason for the overall loss of 
traciality of ip for t <1. 

5. Case = g < t: "t-DEFORMED Free Probability" 

The remainder of the paper considers the case Q = q < t < 1, corresponding to a new single- 
parameter deformation of the full Boltzmann Fock space of free probability [18^119]. Once again, this 
deformation will turn out to be particularly natural. In particular, the statistics of the t-deformed 
semicircular element will be described in an elegant form afforded by the deformed Catalan numbers 
of Carlitz and Riordan [23], the Rogers- Ramanuj an continued fraction, and the t-Airy function of 
Ismail [21]. Moreover, this setting will give rise to a natural counterpart of the result of Wigner (49j . 
as the i-deformed semicircular element will be seen to encode the first-order statistics of correlated 
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Wigner processes. Prior to delving into these properties, we take a moment to show that the von 
Neumann algebra of bounded linear operators on =^o,t(=^) is generated by the creation operators 
{ao,i(/i)}/ie^- 

Lemma 12. For q = < t < 1, the von Neumann algebra #o£ generated by {aot{h)}he^^ is 

Proof. The proof of this fact follows along analogous lines to the full Boltzmann Fock space setting, 
namely the q = 0, t = 1 case. For concreteness, the following sketch adapts the proof of Theorem 
2.2 of In particular, let {ei,e2,...} denote an orthonormal basis for ^ and consider the 

operator P = Z^^i a{ej)*a{ej) G '^o,t- Since 



'0, j 
t""^ei2 (g) . . . (g) ei„, n=i 



a{ej)Cl = and a{ej)ei 



where a{h) := ao^tih), it follows that PCl = and P Cj^ (g) . . . (g e^^ = e^^ . . . (g) ej,^. Therefore, 
P is the projection onto the closure of ©„>i=^^"- Since {^,ri)q^t = whenever ^ G J^®"- and 
r] G J^'^"^ for n ^ m {n,m > 0), P is the projection onto the orthogonal complement of the 
vacuum and Pq := I — P G Wq^i is the projection onto the vacuum. Thus, #o,i contains the 
operator a(eji)* . . . a(ei„)*Pn «(eji) . . . a(ej™), which is a rank-one operator with image spanned by 
(g) . . . (8) ei„ and kernel orthogonal to e^^ (g . . . (g e^^ . It follows that #o,t contains all finite-rank 
operators. Taking closures, #o,t 3 =^(=^o,t(-^)) and the result follows. □ 

5.1. The t-semicircular element. In light of its present interpretation in the context of de- 
formed free probability, the t- semicircular element is a renaming of the (0, t)-Gaussian element 
so,t(/i) G %,t. Specializing accordingly the spectral properties derived in the previous section, the 
t-semicircular element turns out to encode certain familiar objects in combinatorics and number 
theory. 

Lemma 13. The moments of the t- semicircular element so,t(/i) are given by 

r(^NC2{2n) 

where NC2{2n) denotes the lattice of non-crossing pair-partitions and Cn^ are referred to as the 
Carlitz-Riordan t-Catalan (or, rather, q-CatalaT^ num6ers[22l |23], given by the recurrence 



Cli' = Y.^'''ctl&t„ (33) 

k=l 



with C^*^ = 1. 

Proof. The first two equalities are obtained by Lemma [TO] substituting g = 0. Next, set := 
^y&NC2{'in)^^^^^^^^'' ^^"^ '^o^^ ^^^^ '^0 = 1- For TT G NC2{2n), consider the pair (l,/3) G vr and 
note that, since vr is non-crossing, (1) (3 is even, and (2) the remaining pairs either belong to the 
interval / = {2, . . . ,/3 — 1} or to /' = {/3 + 1, . . . ,n}. Thus, given a nesting in vr, either both 
participating pairs belong to I or both belong to /' or one of those pairs is in fact equal to {l,/3}. 
In the latter case, the nesting must be formed by drawing a second pair from /. Summing over all 



vr G NC2{2n) and conditioning on the choice of 2/3 G {1, 2, . . . , n} recovers the recurrence in (33). 



The corresponding argument is illustrated in Figure M □ 



^In keeping with the original objects in [23] , the g-Catalan numbers of CarUtz-Riordan are more frequently given by 
q^^^'^Cn^''\ i.e. with reversed coefflcients compared to the present definition. 
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Chord C 1, 2k) generating k-1 
additional nestings 

0000l)0-000 

V 2k 2n 

k-1 cliords generating n-k chords generating 
m nestings m' nestings 




Figure 9. Recurrence (33): given vr G NC2{2n) with {1,2k) G vr such that the 
chords contained in the interval I = {2, . . . ,2k — 1} generate m nestings and those 
contained in I' = {2k + 1, . . . , n} generate m' nestings, the total number of nestings 
in vr is given hy m + m' + k — 1. 



The first few even moments of so,t(^) are thus given by ipq,t{sq,t{h)'^) = ll^ll^, 'Pq,t{sq,t{h)^) = 
\\h\\j^{l + t), ipq^t{sq,t{h)^) = + 2t + t'^ + t^). The deformed Catalan numbers of the above 

lemma are known to be related to, among other statistics, inversions of Catalan words, Catalan 
permutations, and area below lattice paths |22j . 

Rather than further considering individual moments, specializing to g = the continued frac- 
tion encoding of the {q, t)-Gaussian measure yields a more fascinating object still. Note that for 
notational convenience, and without loss of generality, the remainder of the section considers the 
normalized t-semicirculars canonically written as so,t '■= •so,t(e), for some unit vector e G J^. 

Lemma 14. The moments of the normalized t- semicircular element so,t o,re encoded hy the Rogers- 
Ramanujan continued fraction as 



n>0 1 



1 



t^Z 



t^Z 



1 

The Cauchy transform of the t-semicircular measure fiQ^t associated with so,i(e) is given by 



1 , , _i 5^ra 



-1 



in V: 



2 



E„>o(-i) 



>0\ ^/ (l_t)(l_42)...(l_fn) 

^n{n — 1) 
(l_t)(l_42)...(l_tn) 



Proof. The continued fraction follows immediately from Lemma 10 by letting q = 0. Since so,i(e) 
is bounded, the measure /xo,j is compactly supported. Therefore, the Cauchy transform of /io,t has 
the power series expansion 



jRZ-ri z z 



where M{z) = J2n>o ^0,^(^0 1)^^ (^-S- |3S])- It is well known (e.g. [20]) that the above Rogers- 
Ramanujan continued-fraction can be written as 



„2 



Z]n>o( 1)" {l~t)(l~t'^)...{l-t'^) 

V ( fn(n-l) 

l^n>0\ ^) (i_t)(i_t2),..(i_tn) 



Recalling that the continued fraction encodes M(z) and performing the required change of variable 
yields the desired Cauchy transform. □ 
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In turn, the functions governing the numerator and denominator of the Cauchy-transform of //o,f 
turn out to have a fascinating interpretation as single-parameter deformations of the Airy function, 
a fact discovered by Ismail [21]. Specifically, let the t-Airy function be given by 

Ai(z) = y- A-zT. (34) 

n>0 ^ ^ ^ 



In other words: 
Corollary 1. 



, . . 1 Mi/z) 



r, zAt{l/{zt)) 



„2 



Remark 7. While the function Aq{z) = Y2n>o (i-q)^ (i-q") features at various points in 

Ramanujan's work, the nomenclature is more recent. It was proposed by Ismail [21j upon discovering 
the fact that, analogously to the Airy function in the classical case, the function Aq is involved in 
the large degree Plancherel-Rotach-type asymptotics for the q-polynomials of the Askey-scheme. 

Prior to seeking to further characterize the measure /io,t, it is instructive to take a moment to 
consider its associated orthogonal polynomial sequence. Since the moments of /xo,t are given by 
t-deformed Catalan numbers, the q = subfamily of the {q, t)-Hermite orthogonal polynomials 
of Definition [5] can thus also be viewed as a t-deformed version of the Chebyshev II orthogonal 
polynomials. 

Definition 6. The t-Chebyshev II orthogonal polynomial sequence {Un{z; t)} n>o is determined by 
the following three-term recurrence: 

zUn{z]t) = Un+i{z;t)+e-^Un-i{z;t), 

with 

Uo{z]t) = 1, Ui{z-t) = z. 



Remark 8. The orthogonal polynomial sequence encoded by the generalized Rogers-Ramanujan 
continued fraction was first considered by Al-Salam and Ismail in [27\, and the polynomial Un{z; t) 
is the special case of their polynomial Un{z; a, b) for a = 0,b = t. Similarly to the previous section, 
it again follows from the classical theory [17j that the t-semicircular measure //o,t is the unique 
positive measure on M orthogonalizing the t-Chebyshev II orthogonal polynomial sequence. 

The elegant form of the Cauchy transform of the previous lemma provides means of describing 
this measure via the zeros of the t-Airy function. This was indeed done in [27j and, adapted to the 
present setting, is formulated as follows. 

Lemma 15 (Corollary 4.5 in [27J). Let {zjjjgpj denote the sequence of zeros of the res caled t-Airy 
function At{z/t). The measure /xo,t is a discrete probability measure with atoms at 

with corresponding mass 

2z,A[{zj/ty 

where A[{z) := £At{z). The unique accumulation point of /io,t is the origin. 
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5.2. First-order statistics of the Wigner process. For the purpose of this section, consider 
a Wigner matrix Wn = [wij] G AiNiLoo-{^, ^,^)) to be a self-adjoint random matrix with 
elements {wij}i<i<j<N forming a jointly independent collection of centered random variables with 
unit variance and uniformly bounded moments, i.e. for all i,j E [N], K{wij) = 0, Edwjjp) = 1 
and Edii^jj l"") < Cn G M for all n S N. The asymptotics of the distribution of Wn remain an object 
of extensive study, taking root in the work of Wigner [49j and further evolving over the following 
decades. In particular, considering an expectation functional ipjy : A^Ar(Loo- (K, P)) — )• K given 
by = ;^Tr E, the following result is known as Wigner's Semicircle Law. 

Theorem (e.g. Theorem 2.1.1 and Lemma 2.1.6 in [50]). For Wn a Wigner matrix, the empirical 
spectral measure ofW^ ■= Wn/N converges in probability to the semicircle distribution, with 



hm <^^(M^^n-l) ^ Q 

lim ipNiW^"" 



1 /2n 



Af-s>oo n + 1\ n 



In particular, converges in moments to the familiar semicircular element so,i. It is natural 
to expect that by introducing some correlation into Wigner's framework, some deformation of the 
semicircle law, with an analogous refinement of the Catalan numbers, may be achieved. Indeed, 
this was obtained by Khorunzhy [21] and further developed by Mazza and Piau [2^ in relation 
to the following setup. Let a Wigner process refer to the sequence of {Wj\f,p{k)}k£N of Wigner 
matrices satisfying the following conditions: 

• The moments of W]\f,p{k) are uniformly bounded in k, i.e. for all n G N and all i,j G [N], 
lK{\wij{k)\^) < Cn S M, where Wij{k) corresponds to the {i,jY^ entry of the matrix W]\f^p{k) 
and Cn does not depend on k. 

• The processes {wij{k))keN for i < j form a triangular array of independent processes. 

• Each process {wij{k))keN is p-correlated, i.e. for some \p\ < 1 and any 1 < k < m, 

EK,-(A;)u;„-(m)) = p™-^ (35) 

Note that for p = 0, {Wiyflip)} is a sequence independent, identically distributed Wigner matrices, 
whereas for p = 1, the situation reduces to having copies of the same matrix with W]\f^i{l) = 
WN,i{k) for all keN. Let 



n,N 



N 



N 



N 



The convergence in N of the sequence Bn,N was previously established in [2ll [25] and the corre- 
sponding limits computed explicitly. The surprise is that the limiting moments are, in fact, those 
of the t-semicircular element. 



Proposition 2. For n G N and \p\ < 1, 

'W nA'^) WnM) VFjv,p(2n - 1) 
N N "' N 

WnJI) Wn,p(2) WnA^^) 



lim ifj\f 



lim ifiy 



N 



N 



N 



In particular, for p G (0, 1], 

lim ifj\f 
for t = p2 and \\h\\^ = 



M^7V,p(l) VF^,p(2) ^7V,p(n) 



N 



N 
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■r£NC2(2n) 



2 nest{T) 



(36) 
(37) 



Proof. Expressions (36) and (37) can be recovered from 12A\ I25j (cf. Theorem 1 |25j ) via Lemma 13 
The following self-contained sketch is included for completeness. First note that the general form 
of both expressions follows analogously to Wigner's proof of the Semicircle Law. In particular, 
unrolling the normalized trace tpisf, 

ipNiWNA'^WNA^)---WNAn)) = ^ E nwi„i,{l)wi,,,,{2)...Wi^,i,in)). (38) 

il,...,i„e[iV] 

For a fixed choice of ii, . . . G [A^], the typical argument then proceeds by considering index 
pairs which repeat and viewing the corresponding pattern as a partition of [n]. Since 

the individual elements Wij{k) are centered, any partition containing a singleton, i.e. a block 
formed by a single element, does not contribute to the sum. Next, by counting all the choices of 
indices ii, . . . , z„ S [N] corresponding to a given partition and taking into account the normalization 
factors, it can be shown that only the non-crossing pair partitions contribute in the limit. (Note 
that the counting argument is warranted by the fact that there are uniform bounds, in the time 



variable n, on the higher moments.) This yields (36) and the general form of (37). 



Considering the left-hand side of (38), note that for any index choice ii, . . . ,in G [N] such that the 
repetitions in {11,12}, {i2,h}, ■ ■ ■ , {in, h} encode a pair-partition, the expectation E{wi^^i2 (I)^i2,j3(2) 
• • • '^in,ii ("-)) factors into second moments. Furthermore, recalling that for all i, j, K{wij{k)'Wj^i{m)) = 
pin-k^ one obtains that 

where y = {(ai,6i), . . . , (an,6n)} with 1 = ai < . . . < a„ and m < hi. 

It remains to rewrite the above sum in terms of nestings. But, that 'f is non-crossing implies 
that for any i G [n], p^i~°-^ = pi+2xnest(a,,6,;r)^ where nest(aj, 6^; denotes the number of nestings 
that include the pair {ai,hi), i.e. nest(ai, 6j; 1^) := \{j G [n] | a, < aj < bj < bi}\. Thus, for any 
y G NC2{2n), it follows that . . .p^"""" = p"+2nest(r)^ rpj^-g y-^^^g ^ completes the 

sketch. □ 

Remark 9. Taking a process formed from a copy of the same matrix, viz. WAr.i(n) = Wat a.e. for 
all n G N, yields p = 1 = t and recovers the result of Wigner. In particular, denoting Wj\f = Wn/N, 
the t-Catalan numbers become the (usual) Catalan numbers and 



^' Wn, Wn,i{2) WN,i{n) ^ 
lim 77 77 • • • 77 = lim <fN{WN ") = ¥^0,1(^0,1 



N^oo \ N N N N 



¥00 



i.e. {limAr_j.oo (^jv(I^Af "')}neN is the moment sequence of the semicircular element in free probabil- 
ity. However, for p G (0, 1], it is not obvious that {lim7v->oo Vn(Wn,p{1) Wn,p{2) . . . W]\f^p{n))}neN 
should still form a moment sequence. Yet, this is indeed the case, and the moments turn out to be 
those of a t-semicircular element with t = p^ . A deeper principle underlying this fact is presently 
unclear. Nevertheless, given their ties to an array of fascinating algebraic and combinatorial ob- 
jects, their appearance in relation to limits of Wigner processes, and their fundamental role in 
the generalized non- commutative Central Limit Theorem [26], the {q,t)-Gaussians may harbor an 
additional potential for capturing a broader range of behaviors. 
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